Nota técnica 


Predicciones extremas de lluvia en 24 horas 
en el estado de Zacatecas, México 


* Daniel Francisco Campos-Aranda* e 
Profesor jubilado de la Universidad Autónoma de San Luis Potosí, México 


“Autor de correspondencia 


Resumen 


2014). 
Predicciones extremas de lluvia en 24 horas en el estado de 


Campos-Aranda, D.F (septiembre-octubre, 
Zacatecas, México. Tecnología y Ciencias del Agua, 5(5), 199- 
225. 


Cuando se dimensionan hidrológicamente obras hidráulicas 
grandes o pequeñas pero peligrosas por su ubicación cercana 
a centros de población o de importantes áreas económicas, es 
común utilizar crecientes de diseño asociadas a altos periodos 
de retorno como 1 000 o 10 000 años; incluso en algunas 
ocasiones se decide eliminar el riesgo de ser excedidas 
utilizando la creciente máxima probable. Por lo general, 
la obtención de estas crecientes de diseño no puede ser 
abordada con el método probabilístico, pues no existen datos 
de gasto máximo anual, entonces su estimación se realiza 
con métodos hidrológicos que transforman tormentas de 
diseño en hidrogramas de respuesta. Desafortunadamente, 
los registros pluviográficos son también bastante escasos y 
por ello, se procesan los registros de precipitación máxima 
diaria anual (?MD), que son mucho más comunes y amplios. 
En este estudio se procesaron los 98 registros de PMD con 
más de 25 datos, disponibles en el estado de Zacatecas. Se 
encontró que 17 son no homogéneos y por ello se eliminaron. 
Para el resto, 81 series anuales de máximos, se obtuvieron sus 
predicciones de periodos de retorno elevados variando de 
100 a 10000 años. Se contrastaron las predicciones obtenidas 
con intervalos de recurrencia de 100, 1 000 y 10 000 años, 
con las disponibles como curvas isomáximas publicadas en 
1976, por la extinta Secretaría de Recursos Hidráulicos y se 
encontró que tales mapas todavía son vigentes para el estado 
de Zacatecas. También se aplicó el método estadístico de 
David M. Hershfield para estimar la precipitación máxima 
probable (PMP) puntual en 24 horas. Esta PMP varió de 
196.5 a 507.0 milímetros. Respecto a los cocientes entre la 
PMP y el valor máximo observado de PMD y entre la PMP 
y la predicción de 10 000 años, sus magnitudes medias de 
orden práctico son 4.20 y 2.20, respectivamente. Los análisis 
regionales de frecuencia realizados en dos zonas geográficas 
del estado de Zacatecas, permitieron demostrar su utilidad 
y precisión para estimar predicciones en sitios o localidades 
sin datos, pero ubicados dentro de tal región. 

Palabras clave: precipitación máxima diaria anual, 
precipitación máxima probable, distribución Log-Pearson 
tipo III, distribución GVE, error estándar de ajuste, análisis 
regional de frecuencia. 
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Campos-Aranda, D.F. (September-October, 2014). Extreme 24- 
Hour Rainfall Predictions in the State of Zacatecas, Mexico. Water 
Technology and Sciences (in Spanish), 5(5), 199-225. 


When determining the hydrological dimensions of both large and 
small hydraulic works that are dangerous because of their proximity 
to populations or important economic areas, the use of design floods 
associated with high return periods is common, such as 1 000 or 10 
000 years. At times, the risk to be exceeded is even eliminated by 
using the maximum probable flood. Generally, probabilistic methods 
cannot be used to obtain these design floods since annual maximum 
flow data do not exist. Therefore, it is estimated using hydrological 
methods that convert design storms into response hydrographs. 
Unfortunately, pluviograph data are also scarce and, therefore, 
annual maximum daily rainfall (MDR) records are processed, which 
are much more common and widespread. This study processed 98 
MDR records with over 25 data available for the state of Zacatecas. 
Seventeen were not homogeneous and were therefore eliminated. 
For the remaining 81, annual maximum series predictions were 
obtained for high return periods varying from 100 to 10 000 years. 
Predictions obtained with recurrence intervals of 100, 1 000 and 
10 000 years were compared with those available from maximum 
isocurves published in 1976 by the former Secretary of Hydraulic 
Resources and found that such maps are still valid for the state of 
Zacatecas. The David M. Hershfield statistical method was also 
applied for estimating probable maximum precipitation (PMP) 
punctually with duration of 24 hours. This PMP ranged from 
196.5 to 507.0 millimeters. Regarding the ratios between the PMP 
and the maximum observed value of DMP, and between PMP and 
the prediction of 10 000 years, their practical average magnitudes 
are 4.20 and 2.20, respectively. The regional frequency analyses 
performed in two geographic regions in the state of Zacatecas made it 
possible to demonstrate their usefulness and accuracy in estimating 
predictions for sites or localities in that region that do not contain 
data. 


Keywords: Annual daily maximum  precipitation, probable 
maximum precipitation, Log-Pearson Type II distribution, GEV 
distribution, standard error of fit, regional flood frequency analysis. 
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Introducción 


El análisis de frecuencia de crecientes permite 
estimar los gastos máximos de diseño, que 
son eventos extremos asociados con una cierta 
probabilidad de excedencia, cuyo recíproco 
es el llamado periodo de retorno o intervalo 
promedio en años entre la ocurrencia de un 
evento igual o mayor. Su procedimiento 
consiste en ajustar una función de distribución 
de probabilidades o modelo probabilístico a la 
muestra disponible de gastos máximos anuales 
instantáneos y entonces utilizarla para obtener 
la predicciones buscadas. 

Cuando tales predicciones o crecientes 
de diseño dimensionarán hidrológicamente 
grandes embalses de aprovechamiento o 
control, diques de protección de centrales 
nucleares, o bien buscarán definir la cota 
máxima de inundaciones para ubicar arriba 
de tal nivel, carreteras importantes y plantas 
de tratamiento o de abastecimiento de agua 
potable, los periodos de retorno que se utilizan 
son elevados, variando de 1 000 a 10 000 años, 
e incluso buscando suprimir el riesgo; por ello 
se emplea la creciente máxima probable, la cual 
se estima con base en la precipitación máxima 
probable (Smith, 1993; Gupta, 2008; Shaw et al., 
2011). 

Debido a la escasez de estaciones de aforos 
en los sitios de interés, la estimación de las 
crecientes de diseño se realiza a través del 
llamado método hidrológico, el cual transforma 
una tormenta de diseño en un hidrograma 
de respuesta de la cuenca. El procedimiento 
fundamental de este método es la técnica de 
los hidrogramas unitarios (HU), sean éstos 
identificados o sintéticos. Entre los segundos 
se tienen los HU adimensionales y los 
triangulares, que son su simplificación. 
dada 
pluviógrafos en las cuencas rurales, sean éstas 


Por otra parte, la escasez de 
pequeñas o grandes, el enfoque obligado para 
desarrollar las tormentas de diseño consiste en 
utilizar los registros pluviométricos, que son 
más abundantes y amplios. Los registros de 
precipitación máxima diaria anual (PMD) 


son procesados como variables aleatorias, 
ajustándoles un modelo probabilístico para 
obtener las predicciones buscadas, es decir, se 
desarrolla con ellos un análisis de frecuencia. 

El objetivo de este trabajo consistió en 
analizar y procesar los 134 registros disponibles 
de PMD del estado de Zacatecas, México, para 
estimar las predicciones correspondientes a los 
periodos de retorno de 100, 500, 1 000, 5 000 
y 10 000 años, con base en dos modelos 
probabilísticos que se aplican bajo precepto. 
También se estimó la precipitación máxima 
probable (PMP) puntual en 24 horas de 
duración, por medio del método estadístico 
de David M. Hershfield. Existiendo mapas 
del país de curvas isomáximas de 100, 1 000 y 
10 000 años de periodo de recurrencia, 
publicadas en 1976 por la extinta Secretaría de 
Recursos Hidráulicos (SRH), se contrastaron 
sus valores con las predicciones obtenidas. Se 
encontró que tales mapas todavía son válidos. 
Respecto a las relaciones entre la PMP y el valor 
más grande observado de PMD y aquélla de 
periodo de retorno 10 000 años, se definieron 
los siguientes valores promedio de carácter 
práctico: 4.20 y 2.20, respectivamente. Se 
exponen también los resultados de dos análisis 
regionales de frecuencia, uno para la zona 
geográfica de Juchipila y otro para la Región 
Hidrológica Núm. 37 (El Salado) del estado de 
Zacatecas, México. 


Datos, métodos y resultados 
Registros pluviométricos procesados 


El estado de Zacatecas dispone actualmente 
de 134 estaciones pluviométricas con registro 
de precipitación máxima diaria anual (PMD), 
según información proporcionada por la 
Dirección Local de la Comisión Nacional del 
Agua (Conagua). Tomando en cuenta que tales 
registros serán procesados probabilísticamente 
para obtener predicciones en altos periodos 
de retorno, se decidió utilizar sólo aquellas 
series con más de 25 datos. En el cuadro 1 se 
muestran las características generales de los 
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98 registros con más de 25 años, así como sus 
propiedades estadísticas insesgadas, definidas 
con la media aritmética y los coeficientes de 
variación, asimetría, curtosis y de correlación 
serial de orden 1. El registro más largo es de 
68 años y corresponde a la estación San Pedro 
Piedra Gorda, pero se tienen además 20 series 
con 50 o más años. Previo a procesamiento 
probabilístico y estadístico de cada registro de 
PMD se revisó cada uno para detectar todos los 
valores superiores a 100 milímetros y solicitar 
a la Conagua de Zacatecas su verificación O 
rechazo, así como el dato correcto. 


Pruebas de calidad estadística 


Para que los resultados de un procesamiento 
probabilístico de una serie de PMD conduzca 
a predicciones confiables, tales datos deben 
proceder de un proceso aleatorio estacionario, 
es decir que no haya cambiado con el tiempo. 
Lo anterior implica que la serie de PMD debe 
estar integrada por valores independientes y 
estar libres de componentes determinísticas 
que la vuelvan no homogénea. Al tomar sólo 
los valores anuales de precipitación máxima 
diaria se tiene la seguridad de que ellos 
serán independientes, pero pueden mostrar 
persistencia, tendencia, cambios en la media o 
bien variabilidad excesiva o escasa. 

Para probar la calidad estadística de 
cada registro se aplicaron las siguientes siete 
pruebas: (1) test de Von Neumann es una 
prueba de no aleatoriedad contra componentes 
determinísticas no especificadas (WMO, 1971; 
Buishand, 1982); (2) la prueba de Anderson 
detecta persistencia a través del coeficiente 
de correlación serial de orden 1 (Linsley et 
al., 1988); (3) la prueba de Sneyers también de 
persistencia es más recomendada en registros 
que no son normales (WMO, 1971); (4) el test 
de Kendall detecta tendencia, sobre todo si 
ésta es lineal (Kottegoda, 1980); (5) el test de 
Spearman también de tendencia es más general 
(WMO, 1971); (6) la prueba de Cramer permite 
comparar medias por subperiodos (WMO, 
1971), y por último, (7) la prueba de Bartlett, 


que busca inconsistencia de la dispersión 
(WMO, 1971; Ruiz-Maya, 1977). 


Resultados de las pruebas estadísticas 


En el cuadro 2 se enlistan las 17 estaciones 
pluviométricas que fueron detectadas como 
no homogéneas debido a que muestran las 
componentes determinísticas que se indican. 
Cuando la persistencia O la tendencia sólo 
fue detectada con una de las dos pruebas 
específicas aplicadas, se designó como ligera. 
Donde el test de Von Neumann mostró pérdida 
de aleatoriedad se indicó al registro como no 
aleatorio. 

Además, 12 registros mostraron sólo 
persistencia ligera, detectada con la prueba 
de Sneyers, tales registros fueron procesados 
y son los siguientes: Ameca La Vieja, Boca 
del Tesorero, Camacho, Cedros, García de la 
Cadena, Genaro Codina, Guadalupe, Juan 
Aldama, Nochistlán, Palmillas, San Pedro 
Piedra Gorda y Villa de Cos. Por lo anterior, 
fueron 81 los registros procesados, cuya 
ubicación se muestra en la figura 1. En las 
figuras 2, 3 y 4 se muestran tres registros no 
homogéneos, que no tienen años faltantes. La 
tendencia lineal mostrada en tales figuras fue 
calculada y probada según criterio de Ostle y 
Mensing (1975). 


Predicciones de altos periodos de retorno 


La definición clásica de probabilidad indica que 
es el cociente del número de casos favorables 
al evento entre el número de casos posibles y 
por ello varía de cero a uno; el primer límite 
implica el evento improbable y el segundo 
el evento seguro. Entonces, un evento anual 
extremo que seaigualado o excedido una vez en 
promedio cada 100 años tiene una probabilidad 
de excedencia de 0.01 y su complemento con 
la unidad será la probabilidad de que no 
sea excedido, es decir 0.99, o bien el 99%. Lo 
anterior implica que el evento extremo cuyo 
intervalo medio de recurrencia es de 10 000 
años, su probabilidad de excedencia es 0.0001, 
esto es muy baja, de sólo el 0.01%. 
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Cuadro 2. Estaciones pluviométricas eliminadas por mostrar componentes determinísticas. 


Estación (clave) 


Componente determinística detectada 


1. Adjuntas del Refugio (2) 


Tendencia y variabilidad. 


2. Chichimequillas (12) 


No aleatorio, persistencia, tendencia y cambio en la media. 


3. Cueva Grande (17) 
4. El Platanito (21) 


Persistencia ligera, tendencia y cambio en la media. 


No aleatorio, persistencia y variabilidad. 


5. El Rusio (23) 


Tendencia. 


6. Felipe Ángeles (29) 


Persistencia ligera y tendencia. 


7. Gruñidora (34) 
8. Malpaso (50) 


Persistencia ligera, tendencia y variabilidad. 


Persistencia ligera y tendencia. 


9. Milpillas de Allende (53) 


Persistencia ligera, tendencia y cambio en la media. 


10. Pajaritos de la Sierra (61) 


No aleatorio, persistencia, tendencia y cambio en la media. 


11. Pastorías (64) 


No aleatorio, persistencia y tendencia ligera. 


12. San Juan Capistrano (81) 


Persistencia ligera y tendencia. 


13. Tepetongo (90) 


No aleatorio, persistencia, tendencia y cambio en la media. 


14. Tierra y Libertad (92) 


Tendencia. 


15. Valparaíso (96) 


Tendencia. 


16. Villanueva (102) 


No aleatorio y persistencia. 


17. Zacatecas (103) 


Tendencia y cambio en la media. 


Explicado el concepto de periodo de retorno 
de los eventos anuales, sus predicciones 
respectivas fueron obtenidas con base en dos 
modelos probabilísticos, cuya aplicación en el 
análisis de frecuencia ha sido establecida bajo 
precepto en Estados Unidos y en Inglaterra, 
éstos son respectivamente las distribuciones 
Log-Pearson tipo III (LP3) y la General de 
Valores Extremos (GVE), las cuales tienen tres 
parámetros de ajuste. 

Ambos 
ajustaron con los dos métodos más comunes 


modelos  probabilísticos se 
y consistentes; para la distribución LP3 se 
usó el método de momentos en el dominio 
logarítmico (WRC, 1977) y real (Bobée, 1975); 
en cambio, para el modelo GVE se empleó 
el método de sextiles (Clarke, 1973) y el de 
momentos L (Stedinger et al., 1993; Hosking éz 
Wallis, 1997). 

En el cuadro 3 se indican para ambos 
modelos su menor error estándar de ajuste 
(Kite, 1977) obtenido y las predicciones 
asociadas con los periodos de retorno de 100, 
500, 1 000, 5 000 y 10 000 años. Se han indicado 
con negritas el error estándar de ajuste (EEA) 


del modelo LP3 cuando proviene del método 
de momentos en el dominio real y el de sextiles 
del modelo GVE, pues tales métodos son los 
que en menos ocasiones conducen al menor 
EEA. La enorme similitud que mostraron los 
EEA de ambos modelos probabilísticos, en 
los 81 registros de PMD procesados, definió 
que no era necesario probar y ajustar otra 
distribución, como por ejemplo la Logística 
Generalizada, recientemente establecida bajo 
precepto en Inglaterra (Mansell, 2003; Shaw 
et al., 2011). En la última columna del cuadro 
3 se indica qué modelo probabilístico condujo 
al menor EEA y por lo tanto qué predicciones 
son las correspondientes a cada uno de los 81 
registros procesados. 

Únicamente en 19 de los 81 registros 
procesados, la distribución Log-Pearson tipo II 
condujo a un menor EEA, como se observa en 
la última columna del cuadro 3. En relación con 
los ajustes de la distribución GVE, la mayoría 
resultó tipo 1 o de Gumbel, con parámetro de 
forma (k) cercano a cero. La tipo IT o de Fréchet 
sólo se presentó en cinco estaciones: Corrales, 
El Cazadero, Mezquital del Oro, Sain Alto 
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Durango 


San Luis Potosí 


5e 60* 


(0) 227 
. 
Zacatecas . 76 
64 49* 


25+ 


Aguascalientes 


Los números de las estaciones 
pluviométricas corresponden a 
los nombres indicados en los 
cuadros 3 y 5. 


Figura 1. Localización geográfica de las 81 estaciones pluviométricas procesadas del estado de Zacatecas, México. 


y Sierra Hermosa; en cambio, la tipo Ill o de San Antonio del Ciprés, San Pedro Piedra 
Weibull se obtuvo en 10 estaciones: García de la Gorda, Tayahua, Tecomate, Trancoso y Vicente 
Cadena, Huanusco, Mesillas, Presa El Chique, Guerrero. 
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Precipitación máxima diaria anual (PMD), en milímetros 


Figura 2. Serie cronológica de precipitación máxima diaria anual en la estación pluviométrica Chichimequillas del estado 


de Zacatecas, México. 
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Figura 3. Serie cronológica de precipitación máxima diaria anual en la estación pluviométrica Milpillas de Allende del estado 


de Zacatecas, México. 
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Contraste de predicciones actuales y del medidas en ciertos intervalos fijos para trans- 
estudio de 1976 formarlas a los valores reales observados de 

lluvia de tal duración. Para que la precipitación 
Weiss (1964) estableció los factores medios por diaria corresponda a la de 24 horas de duración 
los que hay que multiplicar las precipitaciones se debe multiplicar por 1.13. Establecida la co- 
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Figura 4. Serie cronológica de precipitación máxima diaria anual en la estación pluviométrica Zacatecas del estado 


de Zacatecas, México. 


rrección anterior, se procedió a encontrar esta- 
ciones pluviométricas comunes al estudio de 
mediados de los años setenta (SRH, 1976) que 
abarcó todo México y del trabajo que se expone 
para Zacatecas; se encontraron 46, las cuales se 
detallan en el cuadro 4 por regiones hidrológi- 
cas. La Región Hidrológica 12 fue subdividida 
en una zona del oeste, la central (Juchipila) y la 
del este. 

En el cuadro 4 se muestran las predicciones 
tomadas de los mapas de curvas de isovalores 
(SRH, 1976) disponibles sólo en los periodos 
de retorno 100, 1 000 y 10 000 años; además se 
tienen las predicciones adoptadas entre los mo- 
delos LP3 y GVE, previamente corregidas por 
el factor de Weiss y por último se exponen los 
cocientes entre los valores citados; esto es, entre 
el valor reciente y el anterior o de los mapas. 

Lógicamente, se observa mayor dispersión 
en los cocientes del periodo de retorno de 10 000 
años que en aquellos asociados al de 100 


años. Los cocientes mínimos se presentan en 
Camacho y Trancoso, y los máximos en El 
Cazadero y Villa de Cos. En términos generales, 
en las regiones hidrológicas de climas áridos o 
semiáridos, es decir la 36 y 37, los cocientes son 
próximos a la unidad y tienen sus valores más 
bajos en la región hidrológica 12 central. Los 
valores medianos de los 46 cocientes muestran 
gran similitud alrededor de un valor de 0.88, 
por lo cual se concluye que los mapas de curvas 
isomáximas (SRH, 1976) son todavía válidos 
para diseño. Tales mapas se presentan en las 
figuras 5, 6 y 7. 


Estimación estadística de la PMP puntual 


El método estadístico de David M. Hershfield 
para estimar la precipitación máxima probable 
(PMP) puede ser consultado en Hershfield 
(1961, 1965), WMO (1973) y Campos (1998). En 
el cuadro 5 se han concentrado los resultados de 
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Cuadro 4. Contraste de predicciones de lluvia en 24 horas con periodos de retorno de 100, 1 000 y 10 000 años, 
en 46 estaciones pluviométricas comunes a los estudios de 1976 y actual. 


Boletín de tormentas Estimación probabilística Cocientes: VR/VA* 
Núm. | Estación pluviométrica 
100 1000 [ 10000 | 100 [1000 [ 10000 | 100 | 1000 | 000 100 1000 | 10000 100 1000 | 10000 
Región Hidrológica núm. 36 
1 Cedros 125 170 225 77.5 108.6 146.8 0.620 | 0.639 0.652 
2 |Juan Aldama 100 125 175 106.3 130.4 150.7 1.063 | 1.043 0.861 
3 San Gil 115 155 190 91.8 118.0 141.9 0.798 | 0.761 0.747 
4 Río Grande 105 150 180 118.2 172.8 239.4 1126 | 1.152 1.330 
5 El Cazadero 100 125 175 120.1 181.3 262.0 1201 | 1.450 1.497 
6 El Arenal 125 170 200 126.7 168.7 212.8 1014 | 0.992 1.064 
7 Sain Alto 115 155 200 119.0 182.7 272.6 1035 | 1.179 1.363 
8 El Sauz 100 125 150 83.2 104.9 125.5 0.832 | 0.839 0.837 
9 Presa Santa Rosa 100 125 150 102.2 142.3 189.5 1022 | 1.138 1.263 
Región Hidrológica núm. 37 
10  |Concepción del Oro 100 150 180 99.2 131.6 164.8 0.992 | 0.877 0.916 
11  |Camacho 125 175 215 67.5 94.6 128.9 0.540 | 0.541 0.573 
12  |Nuevo Mercurio 100 150 190 104.5 161.3 238.0 1.045 | 1.075 1.253 
13 [San Tiburcio 100 125 180 104.2 140.6 179.0 1042 | 1.125 0.994 
14 [San Benito 130 200 250 130.4 179.6 231.9 1.003 | 0.898 0.928 
15 | Agua Nueva 100 125 150 84.9 111.3 138.4 0.849 | 0.890 0.923 
16 [Guadalupe 100 140 175 100.7 125.8 149.2 1.007 | 0.899 0.853 
17 [Cañitas de Felipe P. 100 150 200 104.3 148.8 200.9 1.043 | 0.992 1.005 
18 [Villa de Cos 100 125 150 118.3 151.4 182.7 1183 | 1.211 1.218 
19  |Fresnillo 100 125 150 101.9 128.8 (553 1019 | 1.030 1.035 
20  |Calera de Víctor R. 100 125 150 87.9 108.7 127.7 0.879 | 0.870 0.851 
21  |San Antonio del Ciprés 100 125 150 95.6 113.2 126.9 0.956 | 0.906 0.846 
22  |Villa Hidalgo 125 200 250 113.1 145.1 175.3 0.905 | 0.726 0.701 
23  |Pinos 105 150 190 111.8 142.0 170.7 1.065 | 0.947 0.898 
Región Hidrológica núm. 11 
24  |Sombrerete 130 175 220 94.7 124.6 155.3 0.728 | 0.712 0.706 
25  |[Chalchihuites 150 200 275 111.8 150.6 193.1 0.745 | 0.753 0.702 
Región Hidrológica núm. 12 (oeste) 
26  |Jiménez del Teúl 125 150 200 92.1 124.0 158.2 0.737 | 0.827 0.791 
27  |San José de Llanetes 120 150 200 98.0 147.6 216.7 0.817 | 0.984 1.084 z 
28  |Achimec 125 150 200 102.0 129.6 155.6 0.816 | 0.864 0.778 2 
29  |Monte Escobedo 25 150 200 92.2 108.5 121.1 0.738 | 0.723 0.606 E 
Región Hidrológica núm. 12 (central) 3 
30  |Boca del Tesorero 100 150 160 98.3 127.4 157.2 0.983 | 0.849 0.983 E 
31  |Presa Palomas 100 150 160 101.7 129.0 157.0 1.017 | 0.860 0.981 E 
32  |Tayahua 100 150 160 87.8 102.2 113.1 0.878 | 0.681 0.707 E 
33  |Presa El Chique 100 150 160 85.8 99.7 110.5 0.858 | 0.665 0.691 E 
34  |Tlaltenango 125 150 200 106.9 129.5 149.6 0.855 | 0.863 0.748 E 
35  |Excamé 125 150 200 106.9 127.8 145.3 0.855 | 0.852 0.727 ? 
36  |La Villita 125 150 200 97.0 113.6 129.2 0.776 | 0.757 0.646 s 
37  |Tecomate 125 150 200 88.1 102.6 113.9 0.705 | 0.684 0.570 al 
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Cuadro 4 (continuación). Contraste de predicciones de lluvia en 24 horas con periodos de retorno de 100, 1 000 y 10 000 años, 
en 46 estaciones pluviométricas comunes a los estudios de 1976 y actual. 


38  |Teúl de González O. 125 150 200 113.0 141.9 169.2 0.904 0.946 0.846 
39  |Juchipila 125 150 200 105.7 145.1 192.4 0.846 0.967 0.962 
40  |Nochistlán 125 150 200 113.0 141.3 166.3 0.904 0.942 0.831 
Región Hidrológica núm. 12 (este) 
41  |Trancoso 120 150 170 78.3 82.9 85.0 0.653 0.553 0.500 
42  |Ojo Caliente 125 170 200 100.3 127.7 154.0 0.802 0.751 0.770 
43 San Pedro Piedra G. 100 125 0.553 
44 | Mesillas 100 130 0.860 
45  |Loreto 100 150 1.059 
46  |Villa García 120 150 0.732 
Valor mínimo 100 125 0.500 
Valor máximo 150 200 1.497 
Valor de la media aritmética - - 0.879 
Valor mediano muestral - - - - - - 0.893 0.873 0.848 
Valor mediano poblacional = = = = = = 0.898 | 0.878 0.861 


valor reciente (VR), valor anterior (VA). 
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Figura 5. Curvas isomáximas de periodo de retorno de cien años del estado de Zacatecas (SRH, 1976). 
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Figura 6. Curvas isomáximas de periodo de retorno de mil años del estado de Zacatecas (SRH, 1976). 


la aplicación de este método en los 81 registros 
de PMD. Se observa que las estimaciones de 
PMP muestran gran variabilidad, fluctuando 
de 196.5 mm en Camacho a 507.0 mm en Sierra 
Hermosa; ambas estaciones pertenecientes a la 
Región Hidrológica núm. 37, sin embargo su 
altitud es bastante diferente, con 1 658 msnm 
en Camacho y 2 100 msnm en Sierra Hermosa. 


Cocientes de la PMP puntual 


En las dos últimas columnas del cuadro 5 
se presentan los cocientes de la PMP a la 
máxima precipitación diaria anual observada 
(Po) y a la predicción de PMD de 10 000 


años (P,,). La media aritmética, la moda y la 
mediana poblacional del primer cociente son, 
respectivamente, 4.177, 4.117 y 4.157; en cambio, 
para el segundo cociente son 2.158, 2.044 y 
2.120. La moda y la mediana se obtuvieron 
con base en la distribución Gamma Mixta. 
Entonces un valor de 4.20 veces la Po y otro de 
2.20 por la P,, conducirá a estimaciones rápidas 
y aproximadas de la PMP puntual de 24 horas 
de duración. Para estimaciones más precisas 
se pueden construir curvas de isovalores de 
PMP, o bien realizar un transporte ponderado 
al sitio de interés, empleando las magnitudes 
del cuadro 5 que sean cercanas y pertenezcan a 
la misma zona geográfica. 
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Figura 7. Curvas isomáximas de periodo de retorno de diez mil años del estado de Zacatecas (SRH, 1976). 


Estimación de la PMP para un área y diversas 
duraciones 


Ha sido establecido que los valores puntuales 
de PMP obtenidos al aplicar el método de 
Hershfield a un registro pluviométrico de 
PMD son representativos de un área de 25 km?. 
Entonces, para estimar la PMP correspondiente 
a una cierta área o cuenca (A) mayor de 25 
km”, se debe aplicar un factor de reducción 
por área (FRA), ya que todas las tormentas son 
menores conforme abarcan más extensión. La 
WMO (1973) presentó unas curvas promedio 
de reducción por área, obtenidas a partir de 
análisis P-A-D de grandes tormentas de tipo 
general. La curva relativa a la duración de 24 


horas ha sido representada por el siguiente 
polinomio de 4” grado: 


FRA =1.012021-5.985305-10*- A +1.39148 107% - 4? 
-1.548155-107-4*+6.121556-10% - A* (1) 


en la cual, el área A está en km?. La expresión 
anterior se obtuvo con base en 11 puntos y tiene 
un coeficiente de determinación de 0.99874 y 
un error estándar de la estimación de 0.00257 
(Campos, 1998). 

En cuencas grandes, donde es común tener 
dentro de ella varias estaciones pluviométricas, 
primero se debe obtener el valor ponderado 
de PMP puntual, por ejemplo, por medio del 
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Cuadro 5. Resultados de la aplicación del método estadístico de D. M. Hershfield para estimar la PMP (mm) 


en 24 horas puntual, en las 81 estaciones pluviométricas indicadas del estado de Zacatecas. 


Medias (mm) D. E.' (mm) E F? PMP? Cocientes* 
Núm. [Estación pluviométrica = = 
. am S, un | Km, | Km | (mm) | PMP/Po | PMP/P,, 
1 Achimec 44.668 | 43.896 | 14.494 | 13.566 | 2.846 | 17.565 | 344.7 4.178 | 2.215 
2 Agua Nueva 36.316 | 35.371 | 12.174 | 10.606 | 3.831 | 17.968 | 278.9 3.670 | 2.015 
3 Ameca La Vieja 44.467 | 43.483 | 12.909 | 11.938 | 2.473 | 17.534 | 330.0 4.521 | 2.507 
4 Boca del Tesorero 43.595 | 42.702 | 13.170 | 11.902 | 3.302 | 17.610 | 309.9 3.779 | 1.972 
a Calera de Víctor Rosales 40.998 | 40.378 | 12.095 | 11.369 | 2.781 | 17.736 | 295.0 4.097 | 2.310 
6 Camacho 28.848 | 27.843 8.674 6.740 | 4.622 | 18.327 | 196.5 3.331 | 1.524 
y) Cañitas de Felipe Pescador | 37.476 | 36.433 | 15.357 | 14.185 | 2.719 | 17.903 | 370.4 4.939 | 1.844 
8 Coapas 41.634 | 40.980 | 11.667 | 11.028 | 2.432 | 17.688 | 293.8 4.333 | 2.437 
9 Chalchihuites 43.653 | 42.330 | 16.837 | 14.343 | 4.335 | 17.634 | 360.6 3.451 | 1.867 
10 [Cedros 32.834 | 32.086 9.814 8.784 | 3.235 | 18.124 | 241.7 3.995 | 1.647 
11  |Col. González Ortega 47.770 | 46.153 | 19.925 | 17.573 | 3.406 | 17.428 | 449.1 4.237 | 1.059 
12 [Concepción del Oro 37.300 | 36.304 | 15.627 | 14.213 | 3.292 | 17.929 | 356.2 4.286 2.162 
13 [Corrales 44.027 | 42.525 | 17.048 | 14.937 | 3.319 | 17.592 | 395.7 4.296 1.433 
14 |El Arenal 48.962 | 47.618 | 19.395 | 17.723 | 2.956 | 17.364 | 448.3 4.483 | 2.107 
15 |El Cazadero 42.905 | 41.931 | 17.056 | 15.596 | 3.435 | 17.666 | 386.1 4.043 | 1.473 
16 |El Nigromante 42.786 | 41.493 | 15.628 | 14.247 | 2.633 | 17.627 | 384.6 4.868 | 2.416 
17 |El Romerillo 46.427 | 45.331 | 13.734 | 12.572 | 2.614 | 17.444 | 345.4 4.417 | 2.373 
18 |El Salvador 46.812 | 45.075 | 16.919 | 14.832 | 2.928 | 17.431 | 404.3 4.568 | 1.773 
19 |El Sauz 38.473 | 38.031 | 11.060 | 10.542 | 2.767 | 17.868 | 271.6 4.042 | 2.163 
20 [Espíritu Santo 46.957 | 45.548 | 17.153 | 15.742 | 2.506 | 17.424 | 421.7 4.961 | 1.820 
21  |Excamé 52.759 | 52.309 | 14.089 | 13.713 | 2.165 | 17.195 | 345.4 4.212 | 2.377 
22 [Felipe Ángeles (S) 43.973 | 42.672 | 15.185 | 13.940 | 2.427 | 17.555 | 381.3 4.984 3.454 
23 [Fresnillo 43.900 | 43.102 | 15.519 | 14.506 | 2.971 | 17.609 | 363.2 4.213 2.339 
24 [García de la Cadena 57.700 | 56.754 | 13.680 | 13.019 | 1.962 | 16.888 | 363.2 4.413 | 3.234 
25 [Genaro Codina 43.532 | 42.670 | 10.905 | 10.095 | 2.390 | 17.565 | 288.2 4.314 2.592 
26  |Gral. Guadalupe Victoria 43.460 | 42.277 | 15.828 | 13.854 | 3.842 | 17.633 | 351.3 3.679 1.776 
27 [Guadalupe 43.129 | 41.700 | 15.846 | 13.938 | 3.178 | 17.625 | 374.9 4.359 | 2.513 
28  |Huanusco 48.131 | 47.477 | 12.348 | 11.880 | 1.980 | 17.363 | 323.2 4.552 | 2.946 
29  |Huitzila 61.652 | 60.117 | 16.736 | 15.192 | 2.527 | 16.719 | 416.1 4.224 | 2.379 
30  |Jalpa 50.335 | 49.567 | 10.547 9.696 | 2.695 | 17.257 | 277.2 3.662 | 2.368 
31  |Jerez 40.793 | 40.098 | 11.682 | 10.912 | 2.676 | 17.733 | 289.4 4.176 | 1.952 
32  |Jiménez del Teúl 36.937 | 35.910 | 13.321 | 11.781 | 3.488 | 17.936 | 309.9 4.025 | 1.959 
33  |Juan Aldama 46.824 | 45.818 | 16.170 | 15.303 | 2.234 | 17.436 | 403.4 5.043 | 2.676 
34  [Juchipila 46.547 | 45.795 | 13.206 | 11.977 | 3.707 | 17.488 | 308.1 3.416 | 1.601 
35  |La Florida 44.211 | 43.600 | 11.930 | 11.177 | 2.899 | 17.584 | 291.7 3.838 2.249 
36  |La Villita 53.868 | 53.325 | 11.119 | 10.495 | 2.742 | 17.129 | 282.2 3.437 2.185 
37 [Las Ánimas 42.607 | 41.341 | 11.918 | 10.044 | 3.530 | 17.631 | 286.2 3.727 | 2.442 
38  |Loreto 51.150 | 49.931 | 17.846 | 15.995 | 3.505 | 17.271 | 398.8 3.762 | 1.883 
39  |Los Campos 48.907 | 47.197 | 17.404 | 14.928 | 3.437 | 17.354 | 400.4 4.065 | 1.663 
40 [Luis Moya 51.744 | 50.285 | 17.397 | 15.966 | 2.469 | 17.191 | 429.0 4.783 | 2.591 
41 [Mesillas 42.687 | 41.231 | 14.972 | 12.897 | 3.386 | 17.644 | 351.5 4.140 | 2.724 
42 [Mezquital del Oro 59.041 | 57.346 | 15.218 | 12.767 | 3.521 | 16.860 | 357.4 3.494 | 1.250 
43  |[Momax 46.424 | 45.025 | 14.865 | 13.398 | 2.610 | 17.433 | 370.5 4.631 | 1.462 
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Cuadro 5 (continuación). Resultados de la aplicación del método estadístico de D. M. Hershfield para estimar la PMP (mm) 
en 24 horas puntual, en las 81 estaciones pluviométricas indicadas del estado de Zacatecas. 


44  |Monte Escobedo 49.477 | 48.784 | 11.275 | 10.415 | 2.930 | 17.321 | 280.6 3.538 2.316 
45  |[Moyahua de Estrada 50.884 | 49.913 | 12.193 | 11.117 | 2.706 | 17.229 | 312.7 3.909 | 2.044 
46  [Nochistlán 47.986 | 47.056 | 17.230 | 15.846 | 3.404 | 17.427 | 391.4 3.875 | 2.353 
47 [Nuevo Mercurio 35.263 | 34.162 | 14.864 | 13.406 | 3.121 | 18.019 | 349.9 4.604 | 1.470 
48 [Ojo Caliente 41.632 | 40.706 | 15.817 | 14.548 | 3.182 | 17.719 | 365.1 4.197 | 2.370 
49  |Palmillas 44.489 | 43.046 | 14.208 | 12.308 | 3.165 | 17.542 | 341.5 4.165 | 2.168 
50  |Pinos 48.876 | 48.044 | 15.921 | 14.814 | 3.089 | 17.377 | 370.9 3.954 2.172 
51 [Pino Suárez 48.050 | 46.274 | 17.969 | 15.608 | 3.186 | 17.389 | 418.6 4.360 1.678 
52  |Presa El Chique 44.381 | 43.850 | 11.316 | 10.647 | 2.945 | 17.580 | 278.7 3.706 | 2.522 
53  |Presa Palomas 47.966 | 47.247 | 12.458 | 11.644 | 2.718 | 17.393 | 306.4 3.883 | 1.952 
54 [Presa Santa Rosa 41.210 | 40.590 | 13.816 | 13.033 | 2.947 | 17.740 | 327.8 4.149 | 1.730 
55 [Puerto de San Francisco 43.468 | 42.403 | 11.369 9.279 | 4.591 | 17.616 | 256.5 3.018 | 1.717 
56 [Purísima de Sifuentes 40.257 | 38.181 | 17.169 | 13.447 | 4.322 | 17.787 | 365.6 3.796 | 1.692 
57 [Río Grande 41.811 | 40.472 | 17.270 | 15.447 | 3.206 | 17.703 | 400.2 4.447 1.671 
58  |Sain Alto 44.760 | 42.833 | 15.770 | 12.754 | 3.777 | 17.541 | 354.3 3.893 1.300 
59 [San Andrés 47.319 | 46.014 | 19.406 | 18.015 | 2.608 | 17.432 | 461.4 4.961 | 2.228 
60 [San Antonio del Ciprés 44.905 | 43.867 | 14.117 | 12.803 | 3.002 | 17.538 | 341.0 4.143 | 2.687 
61 [San Benito 46.175 | 44.700 | 20.032 | 18.800 | 2.197 | 17.466 | 493.3 5.736 | 27 
62 [San Gil 35.464 | 34.306 | 14.806 | 13.265 | 3.143 | 18.008 | 350.0 4.605 | 2.466 
63  |San Isidro delos González | 39.327 | 38.306 | 11.685 | 10.267 | 3.282 | 17.800 | 285.4 3.964 | 1.979 
64 [San Jerónimo 40.007 | 39.007 | 12.276 | 11.182 | 2.682 | 17.755 | 310.3 4.497 | 2.613 
65  |San José de Llanetes 38.340 | 37.207 | 12.583 | 11.139 | 3.052 | 17.845 | 308.5 4.333 | 1.423 
66 [San Pedro de la Sierra 44.540 | 43.396 | 13.802 | 12.831 | 2.229 | 17.514 | 356.6 4.953 | 2.606 
67 [San Pedro Piedra Gorda 41.535 | 40.961 | 11.066 | 10.078 | 3.874 | 17.728 | 262.7 3.284 | 2.969 
68 [San Tiburcio 40.135 | 39.028 | 15.419 | 14.067 | 2.913 | 17.775 | 368.3 4.604 | 2.058 
69 [Sierra Hermosa 46.672 | 45.016 | 21.496 | 19.667 | 2.694 | 17.467 | 507.0 5.173 | 1.584 
70  |[Sombrerete 38.650 | 37.237 | 13.602 | 11.579 | 3.417 | 17.835 | 320.9 4.178 | 2.067 
71  |Tayahua 43.939 | 43.074 | 12.314 | 11.027 | 3.448 | 17.592 | 290.9 3.587 | 2.572 
72  [Tecomate 45.962 | 45.339 | 11.333 | 10.550 | 2.954 | 17.497 | 280.0 3.660 | 2.458 
73  |Teúl de González Ortega 51.566 | 50.602 | 15.584 | 14.381 | 2.948 | 17.236 | 366.7 3.943 | 2.168 
74  [Tlachichila 54.680 | 52.883 | 14.556 | 11.700 | 3.839 | 17.055 | 333.4 3.409 | 1.558 
75 [Tlaltenango 53.143 | 52.433 | 13.599 | 12.633 | 3.148 | 17.175 | 325.9 3.535 | 2.178 
E 76 [Trancoso 44.004 | 43.523 | 12.208 | 11.797 | 2.202 | 17.589 | 304.8 4.386 | 3.587 
3 77 [Vicente Guerrero 47.408 | 46.425 | 12.597 | 11.849 | 2.074 | 17.366 | 334.4 4.710 | 3.060 
3 78 [Villa de Cos 49.276 | 48.118 | 17.678 | 16.065 | 3.242 | 17.354 | 400.9 4.001 | 2.194 
E 79 [Villa García 47.523 | 46.857 | 15.871 | 15.277 | 2.268 | 17.428 | 381.8 4.685 | 2.607 
? 80  |Villa González Ortega 40.713 | 39.760 | 13.212 | 12.307 | 2.400 | 17.721 | 336.1 4.850 | 2.483 
E 81 [Villa Hidalgo 46.470 | 45.255 | 17.271 | 15.510 | 3.369 | 17.486 | 391.1 4.011 | PILAS 
E Valor mínimo 28.848 | 27.843 8.674 6.740 | 1.962 | 16.719 | 196.5 3.018 | 1.059 
de Valor máximo 61.652 | 60.117 | 21.496 | 19.667 | 4.622 | 18.327 | 507.0 5.736 | 3.587 
3 Simbología: 
> 1 desviaciones estándar. 
E ? factores de frecuencia (Km, = máximo observado y Km = calculado). 
Ss 3 precipitación máxima probable en 24 horas puntual. 
57 4 


Po = valor máximo observado de precipitación máxima diaria anual. 
P,, = precipitación máxima en 24 horas con periodo de retorno (Tr) de 10 000 años. 
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método de las curvas isoyetas, los polígonos de 
Thiessen o alguna otra técnica más adecuada. 

Para distribuir en otras duraciones la PMP 
en 24 horas puntual o de un área, se puede 
utilizar el criterio propuesto por Schulz 
(1976), el cual consiste en trazar en un papel 
logarítmico una recta paralela a la curva 
envolvente de valores máximos mundiales, en 
la duración de 24 horas y el valor de la PMP. 
La curva envolvente es una recta en el papel 
logarítmico, cuya ecuación es (WMO, 1973; 
Linsley et al., 1988): 

Pin = 421.64: D'* (2) 
en donde P,,y es la precipitación máxima 
mundial en milímetros de duración D en 
horas. Por ejemplo, para obtener la PMP en seis 
horas, siendo la de 24 horas de 513.0 mm, el 
procedimiento es el siguiente: 


513.0 513.0 


= = 113.375 
24" 4.52481 


PMP, =113.375-(6) ” =265.5 mm 


Análisis regional de frecuencia en la zona 
geográfica de Juchipila 


El objetivo fundamental del procesamiento 
regional de la información de PMD disponible 
radica en permitir efectuar predicciones 
en sitios sin información o en aquellos que 
fueron descartados porque sus registros no 
eran homogéneos, ubicados 


dentro de tal región o zona geográfica. En 


pero ambos 


general, para poder realizar tales estimaciones 
en localidades sin datos, se deben cumplir 
ciertos requisitos estadísticos, por ello el 
análisis regional de frecuencia (ARF) primero 
verifica que los registros no sean discordantes 
entre ellos, después prueba la homogeneidad 
hidrológica define 
predicciones adimensionales, que son válidas 


regional y luego 


dentro de tal área, de manera que después con 
sólo escalarlas se obtienen las predicciones 


buscadas. Como los análisis posteriores se 
basan en los momentos L y sus cocientes, se 
recomienda consultar a Stedinger et al. (1993), 
Hosking y Wallis (1997), o Campos (2010), para 
disponer de sus ecuaciones de cálculo. 

Para la zona geográfica del estado de 
Zacatecas que es atravesada por el río Juchipila 
en la Región Hidrológica núm. 12 Parcial 
(río Santiago) se detectaron 14 estaciones 
pluviométricas, para las cuales se les calcularon 
sus cocientes de momentos L. Las magnitudes 
de los cocientes citados se tienen en el cuadro 
6, en el cual además se expone el valor 
obtenido para la Discordancia (D,) de la prueba 
estadística respectiva (Hosking é: Wallis, 1997; 
Campos, 2010). Se observa que únicamente la 
estación García de la Cadena es discordante 
con el resto y por ello se elimina del ARF. 

A continuación se probó la homogeneidad 
hidrológica de la región a través del nuevo Test 
de Langpbein (Fill $: Stedinger, 1995; Campos, 
2012), encontrándose que sólo las estaciones La 
Villita y Nochistlán quedan fuera de las curvas 
de control de la prueba, la primera por arriba 
y la segunda por abajo. Como tal número de 
estaciones dispersas es admisible, se acepta que 
la región analizada es homogénea. En seguida 
se aplicó el método regional de las estaciones- 
años (Garros-Berthet, 1994; Campos, 2008) 
a los 13 registros conjuntados, previamente 
escalados por su media aritmética. El nuevo 
registro de valores adimensionales tiene 544 
elementos. En el cuadro 7 se han concentrado 
las predicciones adimensionales obtenidas con 
los dos modelos probabilísticos ya utilizados y 
menor EEA. 

Ahora, en el cuadro 8 se muestran los 
momentos L (Hosking € Wallis, 1997; Campos, 
2010) de los 13 registros procesados en la 
región de Juchipila, para obtener sus valores 
adimensionales con respecto a l, y luego los 
ponderados con respecto a la amplitud de cada 
registro (1). Habiendo obtenido los momentos 
L adimensionales ponderados (le), se aplica la 
distribución GVE, para obtener la predicciones 
adimensionales con tal método (Campos, 2008); 
los resultados se muestran en el último renglón 
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Cuadro 6. Cocientes de momentos L y discordancias de los registros de PMD en dos regiones geográficas del estado 


de Zacatecas, México. 


Estaciones de la Estaciones de la Región 
región de Juchipila ; h, t, de Hidrológica núm. 37 ¡ 5 t, de 
1. Excamé 0.15277 | 0.12216 | 0.04587 | 0.63 |1. Agua Nueva 0.18368 | 0.13450 | 0.16887 | 0.44 
2. García de la Cadena 0.13455 | 0.01033 | 0.16733 | 3.11 |2. Calera de Víctor R. 0.16385 | 0.15786 | 0.16759 | 0.97 
3. Huanusco 0.14845 | 0.06257 | 0.00533 | 1.28 |3. Camacho 0.16171 | 0.20173 | 0.12884 | 1.27 
4. Huitzila 0.15598 | 0.12714 | 0.06746 | 0.35 |4. Cañitas de Felipe P. 0.22715 | 0.21471 | 0.13823 | 0.85 
5. Jalpa 0.11391 | 0.16296 | 0.21525 | 1.19 |5. Coapas 0.15869 | 0.13002 | 0.13348 | 0.63 
6. Juchipila 0.15428 | 0.20720 | 0.16859 | 0.25 |6. Concepción del O. 0.23536 | 0.15926 | 0.10614 | 1.67 
7. La Villita 0.11611 | 0.11242 | 0.13096 | 0.72 |7. El Salvador 0.19906 | 0.21614 | 0.18622 | 0.80 
8. Los Campos 0.19337 | 0.22281 | 0.21185 | 1.68 |8. Fresnillo 0.19482 | 0.12217 | 0.16558 | 0.47 
9. Mezquital del Oro 0.13810 | 0.27106 | 0.21152 | 1.39 [9. Gral. Guadalupe V. 0.19497 | 0.20525 | 0.20923 | 1.35 
10. Moyahua de Estrada. | 0.13615 | 0.16536 | 0.10380 | 0.41 |10. Guadalupe 0.20361 | 0.06694 | 0.17305 | 1.82 
11. Nochistlán 0.20191 | 0.13691 | 0.09829 | 1.66 |11. Nuevo Mercurio 0.22853 | 0.26452 | 0.14582 | 1.43 
12. Teúl de González O. 0.16853 | 0.15395 | 0.16774 | 0.41 |12. Pinos 0.18295 | 0.15597 | 0.10857 | 0.16 
13. Tlachichila 0.14511 | 0.21936 | 0.20266 | 0.50 [13. San Antonio del C. 0.17838 | 0.10181 | 0.07163 | 0.73 
14. Tlaltenango 0.14496 | 0.13778 | 0.06988 | 0.43 [14. San Benito 0.24664 | 0.18817 | 0.06252 | 1.81 
15. San Jerónimo 0.17574 | 0.11915 | 0.09293 | 0.32 
16. San Tiburcio 0.21435 | 0.18565 | 0.11888 | 0.47 
17. Trancoso 0.15971 | -0.03596 | 0.09517 | 3.05 
18. Villa de Cos 0.20161 | 0.15431 | 0.12490 | 0.10 
19. Villa González O. 0.18662 | 0.13979 | 0.05099 | 1.37 
20. Villa Hidalgo 0.20729 | 0.15366 | 0.15673 | 0.29 
Cuadro 7. Predicciones adimensionales obtenidas con los métodos regionales indicados en la región geográfica 
de Juchipila del estado de Zacatecas, México. 
Método regional, distribución de Periodos de retorno en años 

probabilidades y método de ajuste 100 500 1 000 5 000 10 000 

Estaciones-años. LP3. Mo. dominio real 1.797 2.078 2.196 2.469 2.586 

Estaciones-años. GVE. Sextiles. 1.820 2.101 2.215 2.464 2.565 

Mo. Ponderados. GVE. Momentos L. 1.843 2.170 2.309 2.625 2.760 


del cuadro 7. Se observa que estas últimas 
predicciones adimensionales son las mayores y 
por seguridad son las adoptadas. 

Por último, en el cuadro 9 se muestran 
las predicciones estimadas con el método 
regional, para las siguientes tres estaciones 
pluviométricas: (1) García de la Cadena, la cual 
fue eliminada conla prueba de discordancias; (2) 
Tecomate, que no se incluyó a propósito desde 
el inicio del ARE, y (3) el sitio de Tepechitlán, 
que no cuenta con estación pluviométrica, 
pero que es una ciudad importante. En los tres 


casos, la magnitud media de la PMD se estimó 
con base en los valores observados cercanos. 
También se puede aplicar una ponderación 
en función de la distancia de cada estación 
pluviométrica al sitio bajo análisis. 

Para fines de comparación entre las 
predicciones históricas (H) y las regionales (R) 


se estimó el error relativo (ER) definido como: 


ed _ (PMD 
(PMD,,) 


1) 


histórico 


Tr 


histórico 1 0 0 ( 3 ) 
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Cuadro 8. Ponderación de los momentos L en los 13 registros de precipitación máxima diaria anual de la región geográfica 


de Juchipila del estado de Zacatecas, México. 
Estación n L, L L le, le, le, le, 
1. Excamé 66 52.759 8.060 0.985 0.1528 0.0187 0.0185 0.0023 
2. Huanusco 36 48.131 7.145 0.447 0.1484 0.0093 0.0098 0.0006 
3. Huitzila 25 61.652 9.617 1.223 0.1560 0.0198 0.0072 0.0009 
4. Jalpa 34 50.335 5.734 0.934 0.1139 0.0186 0.0071 0.0012 
5. Juchipila 59 46.547 7.181 1.488 0.1543 0.0320 0.0167 0.0035 
6. La Villita 53 53.868 6.255 0.703 0.1161 0.0131 0.0113 0.0013 
7. Los Campos 30 48.907 9.457 2.107 0.1934 0.0431 0.0107 0.0024 
8. Mezquital del Oro 27 59.011 8.149 2.209 0.1381 0.0374 0.0069 0.0019 
9. Moyahua de Estrada 31 50.884 6.928 1.146 0.1362 0.0225 0.0078 0.0013 
10. Nochistlán 58 47.986 9.689 1.327 0.2019 0.0277 0.0215 0.0030 
11. Teúl de González O. 44 51.566 8.691 1.338 0.1685 0.0259 0.0136 0.0021 
12. Tlachichila 25 54.680 7.935 1.741 0.1451 0.0318 0.0067 0.0015 
13. Tlaltenango 56 53.143 7.704 1.061 0.1450 0.0200 0.0149 0.0021 
Sumas y promedios 544 52.267 - - 0.1493 0.0246 0.1493 0.0241 


Cuadro 9. Contrastes entre las predicciones históricas y las procedentes de los métodos regionales en dos registros 


de la región de Juchipila, Zacatecas, México. 


Periodos de retorno en años 


Estación TP* 

100 500 1000 5 000 10 000 

García de la Cadena H 88.8 94.0 95.7 98.5 99.4 
X,=57.700 mm 159.3 
60.3 

Tecomate 100.8 
X, =45.962 mm 126.9 
25.9 

Tepechitlán X, =53 mm 146.3 


* tipo de predicción: H = histórica; R = regional; ER = error relativo. 


en donde el error relativo se expresa en 
porcentaje y presenta un valor negativo 
cuando la lluvia estimada con predicciones 
adimensionales regionales resultó menor que 
la calculada probabilísticamente o histórica; 
en cambio, cuando conduzca a un valor 
positivo indica que la predicción estimada 
fue superior a la del registro. Para el registro 
discordante de la estación García de la Cadena, 
los ER variaron del 20 al 60% por exceso con 
respecto a las predicciones históricas, que por 
supuesto se consideran erróneas. En cambio, 
en Tecomate, los ER únicamente fluctuaron 


del 9 al 26%, también por exceso, por lo cual 
las predicciones de los métodos regionales se 
consideran bastante aproximadas. 


Análisis regional de frecuencia en la Región 
Hidrológica núm. 37 


El otro ARF realizado incluyó las 20 estaciones 
pluviométricas del estado de Zacatecas que 
pertenecen a la Región Hidrológica núm. 37 (El 
Salado). Se siguió un procedimiento idéntico 
al descrito en el inciso anterior, por ello en el 
cuadro 6 se muestran los cocientes de momentos 
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L y sus respectivas discordancias (D); se 
observa que resultó discordante solamente el 
registro de la estación "Trancoso. La aplicación 
de la versión corregida del Test de Langbein 
encontró que tal región es homogénea, pues 
ninguna de las 19 estaciones pluviométricas 
que la integran quedó fuera de sus curvas de 
control. 

Ahora en el cuadro 10 se tienen las 
predicciones adimensionales obtenidas con el 
método regional de las estaciones-año, según 
los modelos probabilísticos LP3 y GVE de 
menor EEA, con un registro conjunto de 749 
elementos escalados con la media aritmética. 
En el cuadro 11 se presentan los momentos L 
de los 19 registros de tal región y su respectivo 
escalamiento y ponderación. 

El ajuste dela distribución GVEconelmétodo 
de los momentos L ponderados condujo a las 
predicciones adimensionales que se tienen en el 
renglón final del cuadro 10. Tales predicciones 
resultaron las más elevadas y por seguridad de 
las estimaciones regionales son las adoptadas. 
Por lo anterior, en el cuadro 12 se presentan 
las predicciones estimadas regionalmente en 
las tres estaciones pluviométricas de la Región 
Hidrológica núm. 37, que fueron encontradas 
no homogéneas y también la discordante. En 
todas estas estaciones, la magnitud de la media 
de la PMD se estimó con base en los valores 
cercanos observados. 

El valor del parámetro de forma (k) de la 
distribución GVE, calculado en los dos ARF 
realizados, en la zona geográfica de Juchipila 
y en la Región Hidrológica núm. 37, con los 
métodos de las estaciones-años y de momentos 
L ponderados, siempre resultó cercano a cero 
y positivo, entonces define una distribución 
Weibull, por lo cual no se aproxima al 
encontrado como valor máximo mundial 
por Koutsoyiannis (2004) de -0.150, es decir, 
distribución Fréchet, al procesar conjuntamente 
169 registros de PMD de Europa y Estados 
Unidos, con un número de datos variando de 
100 a 154, para un total de 18 065 estaciones- 
años. 


Conclusiones 


En este estudio de la precipitación máxima 
diaria anual (PMD) del estado de Zacatecas, 
México, quedó manifiesta la importancia de las 
pruebas estadísticas específicas que detectan 
componentes determinísticas en las series 
anuales de máximos, ya que de los 98 registros 
analizados con más de 25 datos (ver cuadro 1), 
17 resultaron no homogéneos y por lo tanto no 
susceptibles de procesamiento probabilístico 
para obtener sus predicciones requeridas en el 
dimensionamiento hidrológico. 

Este nuevo estudio probabilístico de la PMD 
del estado de Zacatecas, que utilizó registros o 
series anuales de máximos, que en promedio 
deben tener de 30 a 39 años más, en relación 
con el elaborado por la extinta Secretaría 
de Recursos Hidráulicos (SRH), con datos 
hasta 1974, sólo mostró predicciones mayores 
en ciertas estaciones, como Río Grande, El 
Cazadero, Sain Alto, Nuevo Mercurio, Villa 
de Cos y Loreto. En otras estaciones, sus 
predicciones resultaron casi idénticas, como 
es el caso de El Arenal, Concepción del Oro, 
San Tiburcio, Cañitas de Felipe Pescador, 
Fresnillo, Boca del Tesorero y Presa Palomas. 
Pero en la mayoría (ver cuadro 4), las nuevas 
predicciones fueron sensiblemente menores, 
por ello se concluyó que los mapas de curvas 
isomáximas de la SRH, publicados en 1976, 
siguen vigentes para el estado de Zacatecas, 
México, y se presentan en las figuras 5, 6 y 7. 

Los resultados del método estadístico 
de Hershfield para estimar la precipitación 
máxima probable en 24 horas puntual (PMP), 
concentrados en el cuadro 5, indican gran 
variabilidad en su magnitud desde valores 
cercanos a 300 milímetros en las regiones 
hidrológicas 36 y 37, hasta las magnitudes 
superiores a los 400 milímetros, comunes de la 
Región Hidrológica 12. Respecto a los cocientes, 
en cada estación pluviométrica, entre la PMP 
y la precipitación máxima diaria observada 
(Po) del registro y la asociada con el periodo 
de retorno de 10 000 años (P,,,), se puede decir 
que resultaron bastante estables o similares, ya 
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Cuadro 10. Predicciones adimensionales obtenidas con los métodos regionales indicados en la Región Hidrológica núm. 37 
(El Salado) del estado de Zacatecas, México. 


Método regional, distribución de Periodos de retorno en años 

probabilidades y método de ajuste 100 500 1.000 5 000 10 000 
Estaciones-años. LP3. Mo. dominio logarítmico 2.050 2.430 2.591 2.961 3.120 
Estaciones-años. GVE. Sextiles. 2.100 2.517 2.692 3.090 3.258 
Mo. Ponderados. GVE. Momentos L. 2.108 2.542 2.726 3.148 3.327 


Cuadro 11. Ponderación de los momentos L en los 19 registros de precipitación máxima diaria anual de la Región Hidrológica 
núm. 37 (El Salado) del estado de Zacatecas, México. 


Estación n L, L L le, le, le, le, 
1. Agua Nueva 43 36.316 6.671 0.897 0.1837 0.0247 0.0105 0.0014 
2. Calera de Víctor R. 51 40.998 6.718 1.060 0.1639 0.0259 0.0112 0.0018 
3. Camacho 31 28.848 4.665 0.941 0.1617 0.0326 0.0067 0.0013 
4. Cañitas de Felipe P. 37 37.476 8.512 1.828 0.2271 0.0488 0.0112 0.0024 
5. Coapas 41 41.634 6.607 0.859 0.1587 0.0206 0.0087 0.0011 
6. Concepción del O. 47 37.300 8.779 1.398 0.2354 0.0375 0.0148 0.0024 
7. El Salvador 25 46.812 9.318 2.014 0.1991 0.0430 0.0066 0.0014 
8. Fresnillo 54 43.900 8.553 1.045 0.1948 0.0238 0.0140 0.0017 
9. Gral. Guadalupe V. 45 43.460 8.473 1.739 0.1950 0.0400 0.0117 0.0024 
10. Guadalupe 31 43.129 8.782 0.588 0.2036 0.0136 0.0084 0.0006 
11. Nuevo Mercurio 38 35.263 8.059 2.132 0.2285 0.0605 0.0116 0.0031 
12. Pinos 5 48.876 8.942 1.395 0.1830 0.0285 0.0134 0.0021 
13. San Antonio del C. 37 44.905 8.010 0.816 0.1784 0.0182 0.0088 0.0009 
14. San Benito 28 46.175 11.389 2.143 0.2466 0.0464 0.0092 0.0017 
15. San Jerónimo 30 40.007 7.031 0.838 0.1757 0.0209 0.0070 0.0008 
16. San Tiburcio 37 40.135 8.603 1.597 0.2144 0.0398 0.0106 0.0020 
17. Villa de Cos 45 49.276 9.935 1.533 0.2016 0.0311 0.0121 0.0019 
18. Villa González O. 31 40.713 7.598 1.062 0.1866 0.0261 0.0077 0.0011 
19. Villa Hidalgo 43 46.470 9.633 1.480 0.2073 0.0318 0.0119 0.0018 
Sumas y promedios 749 39.563 = = 0.1971 0.0323 0.1961 0.0319 


Cuadro 12. Estimación de predicciones en las estaciones pluviométricas indicadas de la Región Hidrológica núm. 37 
(El Salado) del estado de Zacatecas, México. 


Periodos de retorno en años 


Estación 
100 500 1 000 5 000 10 000 
Gruñidora (Xx, = 38 mm) 80.1 96.6 103.6 119.6 126.4 
El Rusio (Xx, =43 mm) 90.6 109.3 117.2 135.4 143.1 
Zacatecas (x, =46 mm) 97.0 116.9 125.4 144.8 153.0 
Trancoso (x, = 44 mm) 92.8 111.8 119.9 138.5 146.4 
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que fluctúan poco y sus valores centrales son 
bastante parecidos (ver cuadro 5), pudiéndose 
recomendar 4.20 para el primer cociente y 2.20 
para el segundo. 

En relación con los análisis regionales de 
frecuencia, éstos permiten verificar o rechazar 
las predicciones obtenidas con registros cortos 
de PMD, así como obtener las de registros no 
susceptibles de procesamiento probabilístico 
por ser no homogéneos y, lo más relevante, de 
sitios o localidades donde no hay información 
pluviométrica. Trabajando por subregiones 
hidrológicas o zonas geográficas, lo más 
probable es que las pruebas de homogeneidad 
regional verifiquen la selección de estaciones 
pluviométricas realizada y que las predicciones 
regionales resulten bastante aproximadas, 
como lo mostraron los dos análisis realizados 
en el estado de Zacatecas: uno para la zona del 
río Juchipila, en la Región Hidrológica núm. 12 
Parcial (río Santiago) y el otro para la Región 
Hidrológica núm. 37 (El Salado). 

Por último, se recomienda llevar a cabo este 
tipo de estudios de la PMD en otras regiones 
hidrológicas o estados del país, ahora que 
existen registros amplios y que los métodos 
regionales han demostrado su utilidad, para 
realizar estimaciones o predicciones confiables 
en sitios sin información pluviométrica, las 
cuales son la base de los métodos hidrológicos 
de obtención de las crecientes de diseño, en 
localidades sin datos hidrométricos. 
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